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Abstract. Whether caused by AGN jets, shocks, or mergers, the most definitive evidence for 
heating in cluster cores comes from X-ray spectroscopy. Unfortunately such spectra are essentially 
limited to studying the emission spectrum from the cluster as a whole. However since the same 
underlying emission measure distribution produces both the observed CCD and RGS spectra, X-ray 
imaging can still provide spatial information on the heating process. Using Chandra archival data for 
a sample of 9 clusters, we demonstrate how imaging data can be used to constrain departures from a 
canonical, isobaric cooling flow model as a function of position in a given cluster. The results of this 
analysis are also shown for the deep archival exposure of the Perseus cluster. Such "heating maps" 
can provide constraints on both the location and magnitude of the heating in the cores of clusters. 
When combined with detections and spectral index maps from low-frequency radio observations, 
these maps can be used to distinguish between different models for heating in these objects. 
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INTRODUCTION 

Since the launch of Chandra and XMM-Newton, it has become evident that the X-ray 
emitting gas in the cores of galaxy clusters is being heated and redistributed by powerful 
radio sources. Images of cluster cores over the past few years show a wealth of structure, 
much of which is produced by interactions between rapidly expanding radio lobes and 
the hot intracluster (ICM) medium surrounding them. In particular, large cavities or 
bubbles are now found routinely in the centers of clusters ifxjj _ Based on studies of cavity 
energetics, the global energy budget in these systems seems more than sufficient to 
balance cooling in some objects |2j,l3|]. The specifics of the feedback loop between AGN 
heating and cooling in cluster cores however is still not understood. In particular, the 
details of where this energy deposition occurs and by what physical mechanism energy 
is transferred to the ICM remain open questions. 

In the canonical cooling flow model where gas is assumed to fully cool from ambient 
values down to very low temperatures, strong X-ray emission lines, such as Fe XVII 
at 15.02 A, would be expected. XMM-Newton RGS spectra of cluster cores however 
show a lack of strong X-ray emission lines from gas at temperatures below a few keV 
[0]. Due to instrumental limitations, XMM-Newton RGS spectra are essentially limited 
to studying the emission spectrum from the cluster as a whole. Fortunately, X-ray 
imaging can sample the cluster temperature structure in cluster cores directly, albeit 
at reduced spectral resolution, since the same underlying emission measure distribution 
or dL/dT distribution produces both the observed CCD and RGS spectra. Using simple, 
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FIGURE 1. Left: Integrated value of a determined within the central cooling region for a small 
sample of clusters with deep Chandra imaging data. This value represents a measure of the total 
heating rate in the core and can be compared to similar estimates based on XMM-Newton RGS 
spectra [5]. Right: Values of a as a function of radius in the Perseus cluster determined from 
archival Chandra imaging data. The increase in a with radius is indicative of an increase in heat 
deposition in the ICM jgj. 



parameterized functions for the dL/dT distribution, we show here that one can map the 
departures from the standard isobaric cooling flow model as a function of position in 
a given cluster. These "heating maps" can provide constraints on both the location and 
magnitude of the heating in the cores of clusters. 



EMISSION MEASURE DIAGNOSTICS 

The lack of strong emission lines from gas at temperatures below ~ T a /3, where T 
is the ambient cluster temperature, in XMM-Newton RGS spectra has become a key 
argument for new cooling flow models which include such additional physics as AGN 
heating, thermal conduction, and star formation. In the standard cooling flow model, 
which assumes no heating, the differential luminosity distribution, i.e. the luminosity 
radiated per unit temperature interval, is proportional to the mass deposition rate, M: 

d ± = -—^—M . (1) 
dT 2 jlm p 

Consequently in this scenario, one should observe equal amounts of emission per unit 
temperature interval, including X-ray lines, from gas down to very low temperatures. 
Based on a sample of 13 clusters, [|4|] find a deficit of emission below about ~ 1/3 of the 
ambient temperature relative to the predicted emission from equation [IJ The observed 
distribution appears to be more consistent with an expression of the form, 

dh 5 k . , jT\ a 

-M(a + 1) - (2) 



dT 2 jim„ \ Tq 



where a ~ 1-2 instead of 0, as expected for the standard cooling flow model. In this 
prescription, a is directly proportional to the amount of heating in the cluster. We have 
implemented the simple, parameterized dL/dT distribution model of equation [2] as a 
source model within both the XSPEC and ISIS spectral fitting packages [0, 01 . This 
technique is similar to that of JUl, although they employ a different formalism for the 
dL/dT distribution. 

With this model, we have used archival Chandra imaging data to measure the under- 
lying dL/dT distribution for a small sample of clusters^. Several of these initial objects 
are common to the original sample of [4]. Data were extracted and fit within the entire 
cooling radius for these objects so these results can be compared with the integrated RGS 
spectra directly. The resulting dL/dT distributions are shown plotted in the left panel of 
Figure [TJ The mean value of a for the sample is < a >~ 1.6. These integrated results 
are completely consistent with the RGS results reported by [4j]. We plan to expand this 
analysis to a larger sample of Chandra and XMM-Newton imaging data. 



MAPPING HEATING IN CLUSTER CORES 

Although fits to the integrated spectra make a nice consistency check, our ultimate 
goal is to actually map the departures from a canonical, unheated cooling flow model 
(equivalent to determining the a parameter) as a function of position in a given cluster. 
Figure [2] shows one such map constructed using the deep, 1 Msec Chandra observation 
of the Perseus cluster. It is immediately obvious that a > over the entire cooling 
region. More strikingly, a rises rapidly immediately outside the edges of the radio filled 
cavities. This observed spatial correlation is strong circumstantial evidence that energy 
associated with the radio cavities is being deposited into the ICM at these locations. Due 
to the episodic nature of AGN activity, the ambient, non-zero value of a may represent 
the time-averaged result of several outbursts, while the observed rise in a at the cavity 
boundaries may be due to ongoing heating associated with the most recent outburst. 

Again, although disconnected from the physics of the heating mechanism, such "heat- 
ing maps" can provide constraints on both the location and magnitude of the heating in 
the cores of clusters. As Figure [2] demonstrates, this technique provides a potentially 
powerful diagnostic for determining the underlying dL/dT distribution in clusters and 
by extension the degree and location of heating in the ICM. Using this technique, we 
intend to construct such heating maps for all clusters in the Chandra and XMM-Newton 
archives with sufficient quality data. Ultimately, these maps can provide constraints on 
AGN feedback-based heating models as well as the energy transfer mechanism itself. 
The complex radio bubble and X-ray cavity morphologies observed in cluster cores 
make it clear that spatial information is key to understanding AGN heating in the ICM. 



1 This initial sample includes Perseus, M87, Hydra A, A478, A262, 2A0335+091, A2052, and A4052. 
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FIGURE 2. A map of heating in the core of the Perseus cluster. The map shows the variation 
of the index a, which measures the departure from simple, isobaric cooling, across the inner 
50 kpc. The value of a at each point in the map is determined by fitting a power-law emission 
measure distribution using the 1 Msec of archival Chandra data on Perseus. Regions of rapidly 
changing a represent sites of heat deposition @]. 
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